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I.  INTRODUCTION 


This  report  deals  with  calculations  and  measurements  performed  to  verify  a 
conceptual  1-0  computational  model  of  a  multi -dri ver,  "non-straight"  shock 
tube.  The  interest  lies  in  the  Army's  perceived  need  for  a  large  blast  ther¬ 
mal  simulator  (LB/TS),  (a)  to  reduce  dependence  on  large  scale,  expensive  HE 
nuclear  simulations,  with  attendant  supply,  weather,  and  coordination  prob- 
Tis ,  and  limitation  to  low  yields,  and  (b)  to  have  available  a  facility 
capable  of  testing  through  the  entire  tactical  nuclear  threat  range.  To  guide 
M.S.  development  and  to  assure  meeting  U.S.  requirements  in  a  cost-effective 
manner,  computations  and  small  scale  experiments  are  necessary  to  assess 
design  and  performance  of  possible  LB/TS  configurations. 

Currently  there  exists  a  number  of  smaller  blast-  and  shock-tubes  located 
the  U.t  ritai  ,  Germany,  and  France.  In  particular,  an  operational 
rqe  bias  simulat  r  (LBS)  facility  located  at  Centre  d1  Etude,  Gramat  (CEG), 

i  r 

F.ance,  %c  incorporates  a  number  of  advanced  design  concepts  and  is  of  a 
unique  multi-driver  design.  The  theoretical  modelling  of  this  facility  has 

3 

been  reported  earlier. 
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We  note  the  interesting  calculations  and  experiments  by  H.O.  Amann,  ’  who 
also  treats  a  multi-driver  configuration  for  a  large  shock  tube:  nozzles  are 
attached  to  high  pressure  bottles  of  equal  length.  In  his  calculations  Amann 
lumps  the  nozzle  openings  into  a  single  opening  and  treats  the  gas  flow  one¬ 
dimensional  ly.  This  model  configuration  is  different  from  that  described  here 
and  the  decaying  blast  wave  simulation  is  not  considered. 

In  the  following  sections  we  will  describe  first  the  experimental  efforts 
and  results,  then  the  computational  effort  and  its  results.  A  comparison  will 
then  be  made  of  the  experimental /computational  results. 


J.R.  Crosner  and  J.B.G.  Monzaa ,  "Lccrge  Diameter  High  Performance  Blast  Sim¬ 
ulators."  Proa.  Fifth  Int’l  Symposium  on  Military  Applications  of  Blast 
Simulation ,  Vol.  1,  Stockholm ,  Sweden,  May  23-26,  1977 . 

2 

S.  Gratias  and  J.B.G.  Monzac,  "The  Large-Scale  Nuclear  Blast  Simulator  of 
the  Gramat  Research  Center-Concept ,  Research,  Performance. "  Proc. Seventh 
Int’l.  Symposium  on  Military  Applications  of  Blast  Simulation,  Medicine  Hat, 
Canada,  July  13-17,  1981. 

3 

A. Mark,  "Computational  Design  of  Large  Scale  Blast  Simulators."  AIAA  19th 
Aerospace  Sciences  Meeting.  St.  Louis,  Missouri,  January  12-15,  1981. 

4h.O. Amann,  "Theoretical  and  Experimental  Investigations  on  the  Driving 
Mechanism  of  a  Shock  Wave  Simulator  of  Large  Cross  Section.  "  Ernest  Mach 
Institute  Freiburg ,  W.  Germany,  Report  No.  2/72  May  1972. 

^H.O. Amann,  "Model  Tests  for  a  Large  Diameter  Simulator  Driven  by  Several 
Generators  Filled  with  Compressed  Air."  Proc.  Fifth  Int’l.  Symposium  on 
Military  Applications  of  Blast  Simulation,  Vol.  1,  Stockholm,  May  23-26, 
1977. 
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II.  EXPERIMENTAL  SHOCK  TUBE  AND  PROCEDURES 


A.  LBS  Model  and  Instrumentation 


The  multi-driver  CEG  facility,  as  depicted  in  Figure  la,  is  modeled  by  the 
2-dimensional,  axi -symmetric  shock  tube  shown  schematically  in  Figure  lb.  The 
circular  cross-section  of  the  driver  in  Figure  lb  at  any  axial  location  is  the 
total  cross-sectional  area  of  the  multi-driver  configuration  at  the  correspond¬ 
ing  location.  The  stepped  driver  then  results.  The  2-dimensional  axisymetric 
LBS  shocktube  was  constructed  to  1/37  scale,  and  is  the  physical  representa¬ 
tion  for  the  computational  model. 

A  conical  nozzle  of  16°  half  angle  represents  the  lumping  of  the  seven  6° 
half  angle  nozzles  of  the  CEG  configuration.  Three  piezoelectric  transducers 
in  the  nozzle  wall,  spaced  at  one  inch  (2.54  cm)  intervals  from  the  nozzle 
opening,  monitor  pressure  development  in  the  nozzle.  The  nozzle  opens  into  a 
10  inch  (25.4  cm)  diameter  driven  tube  which  may  be  closed  off  to  the  atmo¬ 
sphere  by  a  removable  annular  plate  fastened  to  the  nozzle  opening.  This  fea¬ 
ture  permitted  duplicating  the  CEG  facility's  removal  of  their  similar  closure 
plate  for  those  shock  overpresures  where  the  rearward-facing  shock  was  expect- 

2 

ed  to  be  swept  downstream  and  out  of  the  nozzle.  Removal  of  the  annulus 
plate  for  these  levels  allows  air  to  be  entrained  into  the  flow  following  the 
shock,  thereby  extending  the  positive  duration. 

The  test  section/driven  tube  was  instrumented  along  the  wall  with  piezo¬ 
electric  gages  at  one  diameter  intervals  beginning  at  one  diameter  downstream 
of  the  nozzle  opening.  Station  numbers  thus  correspond  to  distance  in  diame¬ 
ters  down  the  driven  tube.  Additionally,  a  pitot  probe  (also  a  piezoelectric 
gage)  was  installed  at  Station  7S,  in  radially  1/4  diameter  and  seven  diame¬ 
ters  downstream  of  the  nozzle  and  across  from  Station  7,  the  test  section  loca¬ 
tion  corresponding  to  that  of  the  CEG  facility.  The  driven  tube  was  a  stan- 
dardlO  inch  (25  cm).  Schedule  40  seamless  steel  pipe,  of  18  feet  (5.5  m)  stock 
length.  The  length  was  left  intact  to  observe  the  lengthening  of  test  time/ 
wave  duration  before  the  end-rarefaction  wave's  arrival  back  at  the  test 
station. 

The  shock  tube  driver  was  pressurized  through  a  regulator  and  manually 
operated  solenoid  valve.  The  driver  pressure  at  diaphragm  burst  was  monitored 
by  a  Bytrex-type  gage  through  an  orifice  tube  communicating  through  the  breech 
plug  to  the  driver  volume.  This  gage  was  connected  to  a  peak  reading  meter 
calibrated  to  read  out  the  pressures  directly.  Various  thicknesses  of  Mylar 
and  aluminum  diaphragms  were  used  to  obtain  a  range  of  driver  pressures.  The 
data  from  13  pressure  channels  were  taken  to  an  analog  tape  recorder  according 
to  a  basic  data  acquisition  scheme.  Figure  2,  which  also  shows  the  scheme  for 
data  reduction. 

The  firing  procedure  called  for  pressurizing  the  diaphragm  to  its  natural 
rupture  point.  Some  questions  about  the  "equilibrium"  state  of  the  driver  gas 
were  raised  in  view  of  noise  and  discrepancies  between  prediction  and  experi¬ 
ment.  Subsequent  tests  using  thermocouples  in  the  driver  to  monitor  typical 
temperatures  during  simulated  runs  will  be  presented.  These  will  demonstrate 
the  negligible  temperature  effects  during  pressurizing  as  well  as  an  axial 
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Figure  1.  Shock  Tube  Configurations 
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3.  Detonator-burst  of  Diaphragm.  Additionally,  tests  were  performed  with 
breakage  of  the  diaphragm  by  standard  detonators  (Reynolds  Industries  RP-87- 
Exploding  Bridge  Wire  Detonators).  The  relevant  Shots  are  16-18,  Figure  7. 
These  were  performed  to  determine  the  detonator's  disturbance  to  the  flow. 

They  anticipate  use  of  detonators  to  rupture  the  diaphragms  in  a  future  multi¬ 
driver  shock  tube  model:  respectively,  7a)  detonator  only,  no  flow;  7b)  deton¬ 
ator  with  a  low  level  shock  pressure,  where  detonator  effects  would  be  most  no 
ticeable;  and  7c)  detonator  with  a  high  level  shock  pressure. 

The  result  in  the  no-flow  situation.  Figure  7a,  shows  a  detonation  shock 
of  smal 1 -impul se  moving  through  the  test  section.  For  the  weaker  flow  situa¬ 
tion,  Figure  7b,  the  detonation  spike  merges  with  the  primary  shock  and  has  no 
separate  identity  beyond  Station  5.  Similarly,  for  the  high  level  shock, 
Figure  7c,  the  detonation  spike  is  detected  only  to  Station  6.  In  both  cases 
detonator  perturbations  appear  no  greater  than  the  noise  appearing  on  the 
traces. 

Additionally,  one  may,  in  Figure  8a,  compare  this  shot,  at  Test  Station  7 
(without  annulus  plate)  with  similar  level,  natural -burst  diaphragm  shots, 
without  and  with  annulus  plate,  respectively.  Figure  8b  and  c.  The  set  of 
results  indicate  that  use  of  a  small  detonator  to  rupture  the  diaphragm  has 
negligible  effect  on  test  station  pressure  level  and  waveform. 

4.  Diaphragm  Blockage  of  Nozzle  Opening.  In  the  shots  at  the  highest 
driver  pressures,  it  was  noted  that  the  recorded  shock  pressures  were 

not  increasing  as  expected.  The  suspicion  was  of  "necking  down"  of  the 
required  thicker  diaphragms  by  the  downstream  nozzle  wall.  As  a  check  a  few 
burst  diaphragms  were  examined  for  effective  open  area  relative  to  the 
(design)  nozzle  throat  area.  This  was  accomplished  by  imaging  the  open  area 
in  parallel  light  and  measuring.  The  numbers  appear  in  the  Table  2,  which 
suggests  a  substantial  blockage  with  the  thicker  diaphragms. 


TABLE  2.  BURST  DIAPHRAGM  OPEN  AREA 

%  Open 


Shot 

VP1 

Matl 

Thickness 

Area 

(mil) 

3 

29 

Alum 

20 

89 

12 

136 

Alum 

84 

70 

14 

137 

A1  urn 

84 

67 

15 

179 

41  urn 

104 

63 

The  downstream  nozzle  throat  section  had  originally  been  relieved  to  accom 
odate  20  mil  diaphragms  (1.900"  vs  1.86"D).  Subsequently,  it  was  bored  out  to 
?.20"D  to  permit  use  of  our  thickest  diaphragms.  The  last  shots  21-23  were 
made  with  this  more  open  downstream  throat.  These  runs  resulted  in  the  higher 
shock  pressures  with  increasing  driver  pressures,  as  would  be  expected. 

The  next  sections  will  describe  the  computer  modeling  efforts  and  will 
present  comparisons  between  computations  and  experiment. 

Text  continued  on  page  31 
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Figure  4.  Driver  Gas  Temperature  during  Pressurizing 
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Results  are  shown  for  a  mid-level  pressure.  Figure  4;  other  results  are 
similar.  The  driver-pressure  record  shows  one  or  more  plateaus,  allowing  gas 
to  equilibrate  or  simulating  the  pause  allowing  tape  recorder  to  come  to 
speed,  then  the  final  pressure  rise  to  and  beyond  the  natural  burst  pressure, 
here  65  atm.  or  950  psi .  The  resulting  temperatures  at  the  three  locations 
are  exhibited.  At  the  station  T2  next  to  the  diaphragm  and  farthest  from  the 
air  inlet  an  initial  sharp  rise  due  to  pressurization  is  terminated  by  cooling 
to  the  walls  and  mixing  with  already  present  gas,  with  more  rapid  cooling  when 
pressurization  is  stopped  (no  fresh  heat  addition)  for  the  tape  recorder.  A 
second  rise  follows  when  pressurization  is  resumed,  to  carry  driver  pressure 
to  beyond  the  natural  burst  point,  whereupon  the  experiment  ceases. 

The  maximum  driver  temperature  thus  should  be ~  12  °K  above  ambient,  for 
this  mid-level  shot.  The  gas  temperatures  nearer  the  air  inlet  stations,  T17 
and  T32,  are  seen  to  be  lower  during  the  filling  process.  At  the  highest 

pressure  run,  the  maximum  temperature  was  40  °K  above  ambient,  while  at  dia 

phragm  burst  the  temperature  was  again-  10-12  °K.  Thus  it  is  concluded  that 
temperature  effects  were  an  insignificant  factor  in  experimental /theoretical 
di screpancies. 

2.  Stagnation  Pressures.  It  is  of  interest  to  examine  the  stagnation 
pressure  records.  Station  7s.  The  stagnation  pressure  offers  a  simple  means 
to  detect  density  differences  in  the  flow  not  seen  by  the  wall  pressure  gage. 
As  was  seen  in  Shot  7,  Figure  3,  and  succeeding  shots,  the  cold  driver  gas 
arrives  at  the  test  station  before  observation  time  terminates.  (Heating  of 
the  driver  gas  to  match  interface  conditions  and/or  to  alleviate  cold  gas 
effects  is  looked  at  in  the  theoretical  section;  experimental  results  will  be 
obtained  and  reported  at  a  later  date.) 

Figure  5  shows  results  at  similar  run  conditions  for  the  stagnation  probe 
located  on  axis  when  the  annulus  plate  is  out  (LBS  Model  "A2")-(a),  and  when 
annulus  plate  is  in  (LBS  Model  "Al")-(b).  The  stagnation  probe  is  moved  next 
to  the  wall  in  (c)  and  (d),  respectively,  for  a  similar  comparison.  On  axis, 
the  stagnation  probe  sees  a  somewhat  different,  less  noisy  region  of  cold  gas 
for  the  closed  annulus  plate  "in"-  configuration  than  for  the  "out"  configura¬ 
tion.  This  result  is  reasonable  since  there  is  no  entrained  air  to  mix  with 
the  driver  gas  flow.  The  differences  are  less  obvious  with  the  stagnation 
probe  next  to  the  wall.  It  is  clear  that  stagnation  pressures  recorded  near 
the  wall  do  not  reflect  the  actual  conditions  a  test  target  would  be  exposed 
to  (were  the  air  entrainment  confined  to  the  tube  periphery  as  in  our  1-D 
modelling).  (The  actual  CEG  facility's  configuration,  however,  would  permit  a 
more  uniform  air  entrainment  over  the  cross  section.) 

A  comparison  of  the  dynamic  pressures,  labelled  Station  7s  -  7,  Figure  6, 
derived  from  the  stagnation  probe  positions  on  axis-(a),  midway  to  axis(the 
normal  position)-(b) ,  and  next  to  wall-(c),  illustrates  the  significant  varia¬ 
tion  in  loadings  with  cold  driver  gas  arrival  at  the  test  station. 


Text  continued  page  22 
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It  may  be  observed  that  for  the  highest  level  Shots  21  and  22  the  pres¬ 
sures  measured  are  a  bit  higher  than  would  be  predi^ed,  opposite  to  the  trend 
for  the  lower  level  shots.  With  an  applications  bias  in  mind,  the  pressure 
measured  here  was  the  mean,  gross  level  characteristic  of  the  wave,  which  for 
these  shots  occurred  some  small  time  after  the  incident  shock,  unlike  the 
earlier  shots.  Immediately  behind  the  shock,  f*o w  properties  should  depend 
only  on  the  shock  Mach  number,  as  indeed  they  do  for  Shots  21  and  22  also. 

The  situation  thus  suggests  that  at  the  “Test  Station-7"  one  is  not  observing 
the  model  tube's  fully  developed  flow  behavior  for  these  highest  pressure 
levels. 

In  Figure  3  we  show  pressures  recorded  over  all  the  stations  for  a  medium 
pressure  run.  Shot  7.  It  illustrates  the  experimental  shock  tube  performance 
and  a  number  of  the  wave  processes  occurring  and  propagating  downstream. 
Attempts  to  follow  individual  and  expected  flow  phenomena  down  the  tube  were 
not  overly  successful,  due  to  complexity  and  lack  of  precise  knowledge  of  the 
interactions  occurring  in  the  stepped  driver. 

In  the  traces,  one  typically  sees,  in  the  diverging  section  (Stations  Nl, 
N2)  and  for  the  first  few  stations  (diameters)  down  the  driven  tube,  the  sharp 
rise  in  pressure  with  shock  passage,  then  a  steep  decay  as  the  backward-facing 
shock  and  a  low  pressure  flow  characteristic  of  jet  overexpansion  passes  over 
the  gage.  The  shock  pressure  is  maintained  for  longer  intervals  as  the  shock 
advances  downstream,  ultimately  taking  on  the  desired  decaying  wave  form  at 
the  Test  Station  7  and  following  stations.  A  considerable  amount  of  electri¬ 
cal  and  other  noise  is  present  in  the  traces.  (Gages  were  not  shock-mounted 
in  the  tube  walls.)  Possibly  masked  are  true  flow-induced  fluctuations  in  the 
p/t  histories.  The  sources  and  effects  of  various  other  waves  in  the  flow 
will  be  discussed  more  fully  in  the  computational  sections. 

1.  Driver  Temperatures  During  Pressurizing.  As  mentioned,  experiments 
were  performed  to  determine  temperature  changes  in  the  driver  during  pressur¬ 
izing,  to  see  if  the  changes  were  significant  and  might  cause  departure  from 
ideal  shock  tube  operation.  Three  thermocouples  were  installed  down  the  axis 
of  the  driver  at  positions  labelled  T2,  T17,  and  T32,  (inches  back  from 
diaphragm  station)  to  place  one  near  the  diaphragm,  one  near  mid-driver,  and 
one  near  the  inlet  air  orifice  at  the  end  wall,  respectively.  These  were 
carried  by  a  rod  fastened  to  a  thick  plate  at  the  diaphragm  station  used  to 
seal  off  the  driver,  thus  permitting  a  simulation  of  the  pressurizing 
operation. 
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distribution  in  temperatures  during  the  pressurizing,  which  was  not  entirely 
expected.  Also,  the  stagnation  probe  at  the  test  station  was  moved  radially 
to  observe  effects  of  annular  plate  removal,  and  records  for  these  will  be 
shown. 

B.  Experimental  Results 


The  run  conditions  achieved  with  the  various  diaphragms  used  are  given 
in  Table  1.  A  range  of  shock  pressures  from  3.1  -  31  psi  (21-214  kPa)  was 
obtained  at  the  test  station.  The  computed  shock  pressure  (ratio)  based  on 
measured  shock  velocity  and  shock  tables  is  also  tabulated,  as  a  check  on  gage 
accuracy  and  to  note  any  significant  departure  from  normal  shock  tube  behavior. 


TABLE  1:  FLOW  CONDITIONS 


Measured 


From  Measured  U 
and  Shock  Tablet 


»  >.3 1 


I  LSI  :  LBS 
SHOT.  IB 
STATION.  I 


TWE.  MSEC 


TEST.  LBS  MODEL  A2 
SHOT.  IB 
STATION.  2 
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III.  ONE-DIMENSIONAL  COMPUTATIONAL  MODELING 

This  section  describes  the  one-dimensional  computational  study  performed 
for  modeling  the  two-dimensional  axisymmetric  experimental  shock  tube.  The 
study  involves  computational  simulations  of  related  configurations,  and  re¬ 
sults  of  the  1-D  computational  predictions  for  the  2-D  axisymmetric  experi¬ 
mental  shock  tube.  All  references  to  computations  or  calculations  imply  usage 
of  the  BRL  Quasi-one-dimensional  (Q1D )  Code. 

A.  BRL  Q1D  Code 

One  of  the  basic  research  efforts  at  the  Ballistic  Research  Laboratory 
(BRL)  is  to  computationally  model  shock  tube  processes.  The  BRL  ^1D  Code  is 
an  adiabatic,  inviscid  eulerian  computer  algorithm  adapted  by  BRL  for  this 
purpose.  There  are  many  references  which  present  good  descriptions  of  the 

code's  computer  algorithm. To  summarize,  the  flow  chart  in  Figure  9  illus¬ 
trates  the  basic  computational  cycle  in  the  BRL  Q1D  Code.  As  shown  in  the  flow 
chart,  the  first  step  is  to  lay  out  a  finite  difference  grid  along  the  shock 
tube.  The  computational  grid  may  be  equidi stantly  partitioned  along  the  tube 
length  or  cells  may  be  clustered  about  a  specified  location  utilizing  the 
hyperbolic  function  incorporated  in  the  code.  Thus,  a  large  number  of  grid 
points  may  be  positioned  where  cross-sectional  area  changes  occur.  Typically, 
the  number  of  spacial  grid  points  varies  from  100-1000. 

Then,  the  initial  conditions  are  set  and  normalized.  The  ideal  gas  equa¬ 
tion  of  state,  equation  1,  and  the  Euler  equations,  equation  2,  are  applied  to 
the  shock  tube  where  P  =  Pressure,  y  is  the  ratio  of  specific  heats,  E  =  Total 
Energy/Volume,  p  =  Density,  U  =  Flow  Velocity,  A  =  Tube  Cross-Sectional  Area, 

T  =  Time,  and  X  =  Distance. 


$ 

G.  Coulter  and  others ,  " Experimental  and  Computational  Modeling  of  Rarefac¬ 
tion  Wave  Eliminators  Suitable  for  the  BRL  2.44m  Shock  Tube,"  US  Army 
Technical  Report  ARBRL-TR-02503,  BRL,  ARC,  MD,  June  1983.  (AD  A131  894) 

7  , 

A.  Mark,  "Computations  of  Shock  Wave  Target  Interaction,  "  US  Army 
Technical  Report  ARBRL-TR-02538,  BRL,  APG,  MD,  December  1983.  (AD  B079  28CL) 
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BASIC  COMPUTATIONAL  CYCLE 


SOLUTION 


Figure  9.  Basic  Computational  Cycle 


P  =  (Y-l)  (E  -  HpV2) 


(1) 


9(PA)  +  _A_  (plJA)  =  0 

3t  3x 

3 ( pUA)  +  3 _  (pu2  +  p^A  .  p  |A  =  o,  and 

3t  3x  3x 

3  (EA)  +  3  uA(E  +  p)  =0, 

3t  3x 


(2-a) 

(2-b) 

(2-c) 


New  variables  at  interior  points  are  calculated  from  a  finite  difference 
formulation  of  the  transformed  Euler  equations.  The  Euler  equations  for  con¬ 
servation  of  mass,  momentum,  and  energy  per  unit  volume  are  solved  in  differen¬ 
tial  form  for  the  field  variables  density,  pressure,  temperature,  total  en¬ 
ergy,  and  a  one-dimensional  component  of  flow  velocity  using  finite  differ- 

O 

encing  techniques  attributed  to  Beam  and  Warming. 

New  boundary  values  are  calculated  for  the  open  end  of  the  shock  tube 
using  the  new  variables  at  the  interior  points  and  conditions  suggested  by 
characteristic  theory.  New  boundary  values  are  calculated  for  the  closed  end 
of  the  tube  using  image  points  and  reflection  principles.  Then  the  computa¬ 
tional  cycle  is  repeated  until  the  desired  time  is  reached,  or  until  some 
convergence  criterion  for  a  steady  state  is  satisfied. 

One-dimensional  computational  modeling  has  its  advantages  and  disadvantag¬ 
es.  Advantages  include  1)  inexpensive  to  run,  2)  capability  to  perform  para¬ 
metric  studies  quickly,  and  3)  usefulness  for  obtaining  good  engineering  ap¬ 
proximations  to  a  wide  variety  of  flow  problems.  The  obvious  disadvantage  is 
that  one-dimensional  analysis  is  necessarily  an  approximation.  No  real  flow 
exists  that  is  truly  one-dimensional.  Further,  there  are  many  problems  which 
are  not  adequately  represented  by  one-dimensional  analysis,  for  example,  flow 
from  a  nonaxi symmetric  shock  tube. 

However,  the  axi symmetry  of  the  experimental  shock  tube  and  related  config¬ 
urations  leads  one  to  expect  good  results  from  a  one-dimensional  analysis  of 
their  flow  properties. 


R.  M.  Beam  and  R.  F.  Warming,  ''An  Implicit  Factored  Scheme  for  the 
Compressible  Navier-Stokes  Equations,"  AIAA  J.16,  393-402,  April  1978. 
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B.  Computational  Modeling  of  Related  Configurations 


First,  the  results  of  computationally  modeling  related  configurations  are 
presented.  This  section  is  included  because,  referring  to  Figure  1,  the  1-D 
computational  model  of  the  experimental  shock  tube  is  complicated  with  various 
area  changes.  In  order  to  understand  the  individual  effects  of  these  area 
changes,  and  to  establish  the  credibility  of  the  BRL  Q1D  Code,  similar  area 
changes  are  modeled  separately  and  compared  to  experiments  where  possible. 

The  area  changes  consist  of  A)  a  converging  nozzle,  B)  a  driver  section  with 
complicated  area  changes,  and  C)  a  diverging  nozzle. 

For  the  converging  nozzle,  experimental  and  computational  data  reported  by 
Gottlieb  are  used  to  substantiate  the  accuracy  of  the  BRL  Q1D  Code.  The  compu- 
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tational  algorithm  employed  by  Gottlieb  is  the  Random  Choice  Method  (RCM). 
Unlike  finite  difference  schemes,  such  as  the  Beam  and  Warming  scheme  in  the 
BRL  Q1D  Code,  viscosity  terms  are  incorporated  in  the  governing  equation  of 
motion  to  deal  with  the  presence  of  discontinuities.  In  the  RCM,  shock  waves 
and  contact  surfaces  are  inherently  discontinuous,  occupying  an  unspecified 
position  between  two  grid  points  in  the  flow  field.  To  first  order  accuracy 
the  location  of  such  discontinuities  is  not  exact,  but  their  average  positions 
are.  Furthermore,  numerical  viscosity  effects  and  truncation  error  accumula¬ 
tion  are  not  inherent  in  the  RCM  algorithm.  Thus,  the  RCM  algorithm  is 
considered  a  good  computational  method  for  validation  of  the  BRL  Q1D  Code. 

A  comparison  of  Gottlieb's  experimental  and  RCM-generated  data11  with  the 
BRL  Q1D  Code-generated  data  for  the  converging  nozzle  reveals  excellent  agree¬ 
ment,  Figure  10.  In  Figure  10,  the  large  decays  in  pressure  are  caused  by  the 
reflection  of  the  rarefaction  wave,  created  upon  burst  of  the  diaphragm,  from 
the  closed  end  of  the  tube.  The  rarefaction  wave  travels  into  the  driver  sec- 
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J.  Glimm ,  " Solution  in  the  Large  for  Nonlinear  Hyperbolic  Systems  of 
Equations ,  "  Communications  in  Pure  and  Applied  Mathematics ,  Vol.  18, 
679-715,  1965. 
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T.  Saito  and  I.  I.  Glass.  "Application  of  Random  Choice  Method  to  Problems 
in  Shock  and  Detonation  Wave  Dynamics , "  UTIAS  Report  No.  240 ,  University 
of  Toronto  Institute  for  Aerospace  Studies,  Downsviev,  Ontario,  Canada, 
October  1979. 

11  Private  Communication  from  Dr.  J.  J.  Gottlieb ,  University  of  Toronto 
Institute  for  Aerospace  Studies  to  Dr.  A.  Mark,  Ballistic  Research 
Laboratory,  March  1983. 


Converging  Nozzle  Shock  Tube 


tion,  reflects  from  the  closed  end,  and  interacts  with  the  converging  nozzle 
such  that  the  rarefaction  wave  is  partly  transmitted  and  partly  reflected. 

The  transmitted  rarefaction  wave  overtakes  the  initial  shock,  causing  a  decay 
in  pressure.  The  reflected  rarefaction  wave  cycles  again  in  the  driver  sec¬ 
tion  in  the  same  manner  as  the  initial  rarefaction  wave.  From  Figure  10,  four 
pressure  drops  are  noted  which  correspond  to  four  cycled  interactions  of  the 
initial  rarefaction  wave  with  the  converging  nozzle. 
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In  another  report,  Gottlieb  presents  RCM-generated  results  for  a  shock 
tube  with  area  changes  in  the  driver  section.  Again,  the  BRL  Q1D  Code  data 
agrees  remarkably  well  with  Gottlieb's  results,  Figure  11.  Figure  11  shows 
fluctuations  or  dips  in  the  pressure-time  histories.  These  fluctuations 
result  from  the  interaction  of  the  rarefaction  wave  with  the  parabolic  area 
changes.  As  the  rarefaction  wave  moves  through  the  convergent  part  of  a  para¬ 
bolic  area  change,  a  reflected  rarefaction  wave  is  produced.  As  the  rarefac¬ 
tion  wave  moves  through  the  divergent  part  of  a  parabolic  area  chanqe,a  reflec¬ 
ted  compression  wave  is  produced,  which  follows  the  reflected  rarefaction 
wave.  These  reflected  waves  appear  as  V-shaped  fluctuations,  or  dips  in  the 
pressure-time  histories.  Gottlieb  and  others  note  the  existence  of  head  los¬ 
ses  for  this  and  other  configurations  where  severe  area  contractions  are  pre¬ 
sent  in  the  driver.  Generally,  head  losses  are  largest  for  sudden  area  ex¬ 
pansions  and  contractions.  The  losses  in  sudden  area  changes  are  explained 
in  terms  of  the  vortex  formation  that  occurs  when  flow  separation  takes  place. 
Energy  must  be  supplied  to  maintain  the  vortex  motion  against  the  action  of 
viscous  forces.  Computationally,  the  inclusion  of  head  losses  cause  a  larger 
decay  rate  and  a  lengthening  of  the  positive  phase  duration.  Head  losses 
become  extremely  significant  for  severe  area  contractions  and  high  local  ve¬ 
locities. 

As  a  final  check  on  the  BRL  Q1D  Code's  ability  to  model  area  changes  a 
diverging  nozzle  configuration  is  examined.  Experimental  data  for  a  diverging 
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nozzle  is  obtained  from  a  report  by  Coulter  and  Bertrand.  Once  again,  the 
BRL  QID  Code  does  an  acceptable  job  of  predicting  the  experimental  results, 
Figure  12.  An  interesting  feature  of  the  diverging  nozzle  pressure-time  his¬ 
tory  is  the  spike  at  the  beginning  of  the  waveform.  The  spike  is  attributed 
to  complex  wave  interactions  that  occur  in  the  vicinity  of  the  exit  of  the  di¬ 
verging  nozzle.  An  exact  and  complete  analysis  of  the  spike  formation  process 
is  beyond  the  scope  of  this  study;  therefore,  it  will  be  reported  on  at  a 
later  date. 
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>T.  J.  Coettlieb  and  others,  " Use  of  Perforated  Plates  in  the  Driver  Section 
of  the  PPL  8-ft.  Shook  Tube  to  Produce  Simulated.  Blast  Waves  'dth  Decaying 
Overpressures , "  1st  report,  15  March  1983,  2nd  report  31  March  19a3,  3rd 
and  final  report  7  April  1983.  University  of  Toronto  Institute  for  Aero¬ 
space  studies,  Ontario,  Canada,  M3H  5T 8. 
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George  A.  Coulter ,  "Shock  Pressure  Increase  in  Converaent  Ducts,"  US  Army 
BRL  Memorandum  Report  So.  1625,  BRL,  AFC,  MD,  January''  1965  (AD  463927). 
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Figure  12.  Diverging  Nozzle  Shock  Tube 


To  recapitul ate,  the  BRL  Q1D  Code  results  for  the  individual  area  changes 
compare  favorably  with  experimental  and  computational  counterparts.  In  antici¬ 
pation  of  the  next  section,  which  deals  with  the  combined  effects  of  similar 
area  changes  in  the  experimental  shock  tube,  one  expects  similar  waveform 
characteristics. 

C.  Results  of  BRL  Q1D  Computations  for  2-D  Axisymmetric  Experimental  Shock 
Tube 


The  1-0  computational  equivalent  of  the  axisymmetric  experimental  shock 
tube  is  depicted  in  Figure  lb.  The  1-0  computational  model  preserves  the 
experimental  shock  tube's  cross-sectional  areas  and  lengths.  Thus,  the  com¬ 
putational  model  has  a  stepped  driver,  a  converging-diverging  nozzle  and  a 
long  driven  section. 

Furthermore,  the  initial  conditions  noted  for  the  experimental  runs  were 
used  as  input  to  the  BRL  Q10  Code,  Table  1.  Using  the  1-0  computational  model 
and  the  initial  conditions  as  input,  the  BRL  Q1D  Code  generates  pressure-time 
histories.  The  pressure-time  histories  become  more  familiar  through  an  examin¬ 
ation  of  the  general  characteri sti cs  of  the  computational  waveforms  and  para¬ 
metric  studies  of  the  pressure-time  histories. 

The  general  characteristics  of  the  pressure-time  histories  at  the  test 
station  are  a  sharp  increase  in  pressure,  the  shock  arrival,  and  a  subsequent 
decay.  The  decay  typically  shows  superimposed  fluctuations.  The  fluctuations 
result  when  the  rarefaction  wave,  produced  upon  breaking  the  diaphragm,  propa¬ 
gates  from  the  diaphragm  location  into  the  driver.  As  the  rarefaction  wave 
moves  through  convergent  area  changes  (stepped  driver),  reflected  rarefaction 
waves  are  produced  which  move  in  the  opposite  direction,  following  the  blast 
wave  front.  These  reflected  rarefaction  waves  show  up  as  small  pressure  drops 
in  the  pressure-time  histories.  Also,  the  waveform  shows  a  significant  pres¬ 
sure  drop  within  the  decaying  signature.  This  large  gradient  is  caused  by  the 
reflection  of  the  rarefaction  wave  from  the  closed  end  of  the  shock  tube.  The 
computational  and  experimental  pressure-time  histories  for  lower  driver  pres¬ 
sures  (Figures  16a  through  d)  show  spikes  at  the  beginning  of  the  decay.  As 
stated  earlier,  this  effect  is  directly  related  to  the  presence  of  the  diverg¬ 
ing  nozzle.  At  the  highest  driver  pressure,  the  waveform  (figure  16c)  shows 
an  abrupt  pressure  drop  which  is  caused  by  the  arrival  of  a  backward  facing 
shock  at  the  test  station. 

The  computational  parametric  studies  performed  for  the  shock  tube  involve: 
A)  temperature  effect  of  pressurization,  B)  heated  versus  unheated  driver,  and 
C)  diaphragm  blockages. 

As  discussed  in  the  experimental  section,  temperature  changes  in  the  driv¬ 
er  are  recorded  during  pressurizing  for  simulated  shots  over  a  range  of  driv¬ 
er  pressure.  One  temperature  distribution  recorded  at  diaphragm  "burst"  is 
input  to  the  BRL  Q10  Code  to  see  if  any  changes  occur  in  the  computationally 
generated  pressure-time  histories.  No  significant  changes  are  caused  by  the 
temperature  distribution.  Therefore,  the  temperature  effect,  caused  by  pres¬ 
surization,  is  concluded  to  be  negligible  for  such  shots  and  not  a  source  of 
discrepancies  between  the  experimental  and  computational  records.  The  tem¬ 
perature  effect  would  become  more  significant  at  higher  pressures;  however,  the 
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Shots  5  and  6,  Figures  16d  and  e  ,  represent  medium  pressure  shots;  driver/ 
driven  pressure  ratios  are  47  and  65.  A  common  characteristic  of  these  shots 
is  the  removal  of  the  annulus  plate  from  the  experimental  shock  tube  to  allow 
mass  entrainment.  Mass  entrainment,  particularly  with  increasing  driver  pres¬ 
sures  sustains  the  positive  phase  pressure  level.  Computationally,  the  mass 
entrainment  is  not  accounted  for;  therefore  in  Figures  16a  and  b,  the  pressure 
level  the  BRL  Q1D  Code  predicts  for  the  flow  behind  the  shock  wave  is  lower 
than  the  levels  recorded  experimentally.  However,  the  comparison  of  trends  in 
the  traces  due  to  rarefaction  waves  is  in  satisfactory  agreement.  Also,  at 
earlier  times  the  pressure-time  histories  agree  extremely  well,  except  for 
small  discrepancies.  The  small  discrepancies  are  attributed  to  experimental 
phenomena,  such  as  Mach  reflections  and  possibly  some  spike  formation  in  Shot 
5. 


Shot  7  and  12,  Figures  16f  and  g,  represent  high  pressure  shots;  driver/ 
driven  pressure  ratios  are  99  and  135.  Shot  7  is,  performed  with  the  annulus 
out,  thus  entraining  mass,  while  Shot  12  is  performed  with  the  annulus  in. 

The  discussion  concerning  mass  entrainment  for  Shots  5  and  6  also  applies  to 
Shot  7.  Upon  comparing  the  experimental  and  computational  pressure-time  his¬ 
tories  for  the  two  shots,  one  notes  the  same  characteristics  in  both  cases. 

The  computationally  predicted  pressure  levels  at  early  times  are  greater  than 
the  experimental  levels.  This  implies  that  head  losses,  discussed  earlier, 
and  other  types  of  losses  (diaphragm,  friction,  heat  transfer,  etc.)  are  sig¬ 
nificant  experimentally  and  should  be  accounted  for  computationally.  Also,  as 
described  in  the  experimental  section,  diaphragm  blockage  effects  are  signifi¬ 
cant  for  these  high  pressure  level  shock  tube  runs.  For  Shot  7,  the  differ¬ 
ence  in  pressure  levels  is  not  too  great  and  the  trends  due  to  rarefaction 
waves  agree.  However,  for  Shot  12,  the  pressure  levels  are  significantly 
different  and  the  decay  trends  disagree.  The  comparison  for  Shot  12  is  un¬ 
satisfactory,  but  by  computational ly  modeling  losses,  currently  under  devel¬ 
opment,  it  is  believed  the  comparison  could  be  improved.  The  computational 
losses  should  lower  the  early  pressure  levels  and  lengthen  the  positive  phase 
duration,  thus  providing  a  better  match  to  the  experimental  data. 

Shot  23,  Figure  16h  is  a  reiteration  of  Shot  12.  However,  instead  of 
computationally  modeling  the  diaphragm  blockage  from  estimates  (Shot  12),  the 
diaphragm  blockage  is  alleviated  as  described  in  the  experimental  section.  A 
comparison  of  Shot  12  and  Shot  23  indicates  a  much  better  experimental/computa¬ 
tional  match  for  Shot  23.  This  implies  the  diaphragm  blockage  is  probably  not 
accurately  modeled  and/or  experimentally  determined  for  Shot  12.  Furthermore, 
Shot  23  is  performed  with  the  annulus  out,  thus  entraining  mass.  This  effect 
is  evident  in  the  longer  positive  phase  duration. 

Finally,  Shot  21  represents  the  highest  pressure  level  shock  tube  run 
(PRAT  =  202).  Again,  as  seen  in  the  waveform  match  before  5  milliseconds, 
Figure  1 6 i ,  losses  are  significant  experimentally  and  are  not  accounted  for 
computationally.  The  arrival  of  the  backward-facing  shock  is  captured  within 
one  half  millisecond.  The  discrepancy  after  5  milliseconds  between  the  wave¬ 
forms  is  attributed  to  mass  entrainment  which  is  not  computationally  modeled. 

In  Fig.  17  we  exhibit  the  overall  performance  of  the  LBS  model  and  its 
comparison  with  the  quasi-ID  calculations.  Other  available  data  from  the 
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Figure  15.  Diaphragm  Blockage  Effects 


Figure  14.  Heated  Driver  Effects  -  Dynamic  Pressure-Time  History 


Figure  13.  Heated  Driver  Effects  -  Pressure-Time  History 


experimental  simulations  showed  a  negligible  temperature  rise,  thus  a  negli¬ 
gible  temperature  effect  for  our  operation. 

The  next  effect  studied  is  heating  of  the  driver  gas.  The  driver  gas  is 
heated  computationally  to  a  temperature  ratio  of  1.53  and  2.0  for  the  initial 
conditions  of  shot  23.  The  resultant  P-t  histories  are  compared  to  the  ori¬ 
ginal  run  (temperature  ratio  of  1.0).  As  seen  in  Figure  13,  heating  of  the 
driver  gas  causes  the  shock  overpressure  to  increase  and  the  decay  rate  to 
increase.  Also,  as  shown  by  Figure  14,  heating  of  the  driver  gas  causes  the 
contact  discontinuity  to  decrease. 

Finally,  diaphragm  blockage  effects  are  studied.  Experimentally,  when  the 
diaphragm  breaks,  its  petals  fold  back  onto  the  shock  tube,  thus  neckinq  down 
the  throat  area.  As  described  in  the  experimental  section,  estimates  of  the 
open  throat  area  are  determined.  Computationally,  the  diaphragm  blockage  is 
modeled  as  a  parabolic  area  contraction.  The  open  throat  diameter  determined 
from  the  measurements  is  input  as  the  throat  diameter  of  the  parabolic  area 
contraction.  Figure  15  presents  pressure-time  histories  with  various  throat 
diameters  for  shot  12.  The  throat  diameters  modeled  are  representative  of  no 
diaphragm  blockage,  the  diaphragm  blockage  approximately  determined  from  the 
measurements,  and  the  diaphragm  blockage  necessary  to  produce  a  shock  over¬ 
pressure  match.  With  no  blockage  (diameter  of  nozzle  throat  =  1.86  inches) 
the  experimental  and  computational  P-t  histories  are  very  different.  The 
match  becomes  better  when  the  blockage  effect,  as  approximately  determined 
from  experimental  measurement,  is  modeled  (diameter  of  nozzle  throat  =  1.66 
inches).  However,  the  best  match  between  the  experimental  and  computational 
P-t  histories  is  for  slightly  more  blockage  than  was  indicated  by  measurement. 

An  11.9%  reduction  (diameter  of  nozzle  throat=  1.46  inches)  in  the  measurement 
diameter  provided  a  good  computational/experimental  match.  The  figure  clearly 
shows  that  diaphragm  blockage  effects  are  important  and  should  be  carefully 
determined  and  modeled. 

D.  Comparison  of  1-D  Computational  Predictions  to  Experimental  Results 

A  comparison  of  the  computationally  generated  P-t  histories  and  the  experi¬ 
mental  P-t  histories  reveals  satisfactory  agreement. 

Shot  20,  Figure  16a,  represents  the  lowest  pressure  shock  tube  run  performed; 
the  <Jri ver/driven  pressure  ratio  (PRAT)  is  6.9.  The  computational /experimental 
match  is  excellent,  even  for  the  spike  at  the  beginning  of  the  waveform, 
which  implies  the  spike  is  basically  caused  by  a  one-dimensional  phenomena. 


Shots  2  and  3,  Figures  16b  and  c,  represent  low  pressure  shots;  driver/ 
driven  pressure  ratios  were  16  and  29  respectively.  Agreement  after  approx¬ 
imately  nine  milliseconds  is  excellent.  Before  nine  millisecond,  the  P-t 
histories  show  a  spike,  that  the  BRL  Q1D  Code  predicts  as  smaller  in  magnitude 
and  duration  than  is  recorded  experimentally.  Since  the  code  can  only  predict 
1-D  phenomena,  it  seems  reasonable  that  there  are  2-D  phenomena,  possibly  Mach 
reflections,  jetting,  etc.  that  occur  during  the  experiment  in  the  diverging 
nozzle  area  to  enhance  the  spike  formation. 
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full-scale  CEG  facility  ’  are  included  for  perspective.  Reasonable  to  very 
good  correlation  is  noted  between  all  the  data  sets  at  the  lower  shock  levels. 
At  the  highest  levels  the  LBS  model  data  is  trending  below  the  Q1D  calcula¬ 
tions.  However,  as  noted  in  the  earlier  sections,  the  calculations  consider 
only  an  ideal  gas  flow;  possible  losses  from  2-D  and  non-ideal  effects,  eg. 
viscosity,  vortices,  etc.,  plus  imperfect  tube  operation,  are  not  accounted 
for.  Within  such  limitations,  one  concludes  that  the  ID  code  predicts  the  LBS 
model's  overall  performance  quite  creditably. 


IV.  SUMMARY  AND  CONCLUSIONS 

An  attempt  was  made  herein  to  verify  predictions  of  the  BRL  Q1D  Code 
against  the  flow  from  a  complicated,  non-straight  shock  tube  configuration. 
With  proper  accounting  for  real  world  behavior,  one  observes  that  the  computa¬ 
tions  and  experiment  agree  reasonably  well  over  the  full  range  of  shock  pres¬ 
sures  with  respect  to  pressure  levels,  pressure  histories,  appearance  of 
spikes  and  fluctuations  induced  by  tube  geometry,  temperature  effects  on  dri¬ 
ver  due  to  tube  operation  procedure,  and  effects  of  throat  blockage  due  to 
diaphragm  opening. 

Thus  the  1-D  code  provides  good  performance  predictions  and  correct  trends 
for  even  very  complicated  flow  geometries,  with  efficiency  and  low  cost. 
Further  computations  utilizing  2-D  codes  and  accounting  for  2-dimensional  and 
real  gas  effects  should  provide  even  better  agreement  but  at  greater  costs  in 
computing  time  and  effort. 
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